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Abstract. We discuss the coupler system of two nonlinear oscillators excited by an 
external coherent field prepared in a maximally entangled state (Bell-like state). We 
show that as a result of the coupler interaction of the system with external broadband 
squeezed vacuum bath, entanglement decay dynamics can be considerably affected. 
Besides the phenomena of sudden entanglement death and its rebirth, a shortening (or 
lengthening) of the total disentanglement time td can be observed, depending on the 
squeezing parameters. Moreover, on the example of one of the reborn entanglement 
cases it is shown that by changing the values of these parameters the maximal values 
of the negativity for the 3 ® 2 system discussed can be tailored. 
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1. Introduction 

Modern quantum physics still takes a very lively interest in quantum entanglement. 
The entanglement problems are exciting not only from the cognitive point of view, but 
also as a source of future applications in the field of quantum information theory. In 
real situations the system which is able to generate entangled states always interacts 
with an external environment, and this interaction leads to the losses of entanglement 
in the systems considered (for instance, for discussion concerning decoherence processes 
for a system of two qubits in a lossy cavity see [TJ and the references quoted therein). In 
particular, if we deal with the cases when the system interacts with a zero-temperature 
reservoir, one can observe asymptotic entanglement decay (if there is no interaction 
between the qubits which produce an entangled state) — meaning that entanglement 
decays to zero for time — > oo. In contradiction, interaction with a thermal reservoir 
is able to shorten significantly that time and entanglement is destroyed in finite time 
(such a behaviour is referred to as the sudden entanglement death [21 El 0]). It is also 
known that when the entangled system interacts with squeezed vacuum bath, sudden 
death of entanglement can be observed [5]. It should be mentioned that under certain 
conditions it is possible to rebuild entanglement in the system considered after its 
sudden death (this is called the sudden entanglement birth [61 [8]). Moreover, such 
"sudden" phenomena can be observed not only for entangling optical systems but also 
for other models, for instance dealing with spins. Thus, quite recently Wang et al [9] 
have discussed sudden spin-squeezing death effect. 

Moreover, to avoid decoherence processes several schemes were proposed. For 
instance, some of them use decoherence free subspaces (DFS) method. In particular, 
one can show that the systems interacting with squeezed environment states taken from 
DFS will never show sudden death phenomena. Contrary, when the interaction involves 
states that are taken not from DFS the system will loose entanglement in a sudden way 
(with possible birth events) [TO] . 

In the present paper we will focus on the decay of initially maximally entangled 
state generated via the coherently excited nonlinear coupler (two nonlinear, Kerr- 
like oscillators interacting nonlinearly with each other) system. Such couplers can 
exhibit numerous interesting physical phenomena not only in the domain of quantum 
information theory but also in other fields of quantum optics (for the thorough review 
of the problem concerning nonlinear couplers see for instance review paper [UJ and the 
references quoted therein, whereas for problems of nonlinear quantum-optical oscillators 
in general see [12J). We shall show that we can manipulate the character of the 
entanglement decay to some extent via the phase of squeezed reservoir in which the 
system is immersed. 
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2. The model and its solutions in squeezed reservoir 

The system we deal with is composed of two nonlinear oscillators labelled by a and b, 
and characterised by the nonlinearity constants Xa and Xb, respectively The oscillators 
are described by the Hamiltonian expressed by the bosonic creation and annihilation 
operators in the following form H n i = ^■{a i ) 2 a 2 + ^(b^) 2 b 2 , where the parameters Xa 
and Xb can be identified as third order susceptibilities of the optical nonlinear media 
- we can refer to our oscillators as "Kerr-like" ones. Such nonlinear systems can be a 
source of not only entangled states but also other relevant states commonly discussed 
in the papers devoted to the quantum optics problems. For instance, Miranowicz at 
al. [13j have shown that Kerr-like oscillator can be a source of discrete superpositions 
of arbitrary number of coherent states (so-called Schrddinger cat-like or kitten states). 
The oscillators interact with each other by an nonlinear coupling of the strength e - the 
Hamiltonian corresponding to this interaction can be written as: Hi = e{a)) 2 b 2 + h.c. 
Additionally one of the oscillators (a) is externally driven by a coherent field. The 
intensity of this linear interaction with the field is equal to a and the Hamiltonian 
corresponding to this process is H e = aa) + h.c. 

While the system is not subjected to the interactions with any kind of external 
environment it can be treated as a quantum scissors device under some assumptions. 
The number of two-mode states engaged in the dynamics can be substantially truncated 
as both internal (characterised via e) and external (described by a) interactions are small 
if compared with nonlinearity parameters x- Therefore, as explained in [H], the number 
of two-mode states effectively engaged in the system's dynamics is drastically reduced. 
The fidelity between such truncated state 1^)^* = c 2|0) a |2) fc + Ci 2 |l) a |2){, + c 2 o|2) a |0)b 

and a two-mode state = J2 c n,m(t) \ n ) a \ m )b * s a l mos t equal to unity with the 

n,m=0 

accuracy ~ 10~ 4 . We show this fact in Fig.l where the fidelity is plotted for e = 0.1, 
a = e and x = 25, and the numerical calculations were performed in the n-photon 
Fock basis where 100 two-mode states were involved. This result enables us to treat our 
system as a nonlinear quantum scissors device. In fact, such behaviour originates from 
the resonances between the energies of the levels generated by the Hamiltonian H n i and 
the couplings discussed in the model. For further informations concerning linear and 
nonlinear quantum scissors devices see for example [151 HI El EES US] and the references 
quoted therein or the review papers [201 EI] • 

As mentioned above, the system can beave as quantum scissors device described 
in details in [H]. This system is able to form Bell-like states that are examples of 
maximally entangled states (MES). It should be mentioned that for the cases when 
interactions with external field are absent, we deal practically with qubit-qubit system, 
whereas when this interaction is present our system should be treated as qutrit-qubit 
one. In consequence, thanks to the two interactions included in our model, we can 
generate the Bell-like states during the system evolution. They can be expressed in the 
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following form: 

\Bi) =^(|2>JO) 6 + i|0)j2) 6 ), (1) 
W =^(|2>JO> k -<|0> 8 |2> 6 ), (2) 

=^(|2>J0> 6 + |1}J2} 6 ). (3) 

Unfortunately, if we include the interaction of the system with external reservoir one 
can expect that such generation process becomes imperfect and we lose the entanglement 
in our system. It is so because the dynamics of the system cannot be restricted to the 
two-mode states which are engaged in \^) cu t any more. The interactions with reservoir 
cause that the other two-mode states (the states with lower energy) can be populated 
and therefore, the system as a whole cannot be treated as a qutrit-qubit one. Moreover, 
due to the interactions with external field (a) also the states with higher energy are 
populated, but to a lower extent, as this interaction is small (in comparison to x value). 
The weakness of this interaction is a necessary condition for the application of the 
nonlinear scissors mechanism. In further considerations we shall concentrate on the 
initially populated subspace defined by the Bell-like state and trace the entanglement 
within this subspace. In such a way we are not interested in the entanglement defined 
for the whole system that can be produced during the interaction time, but in the 
previously occupied 3 <g) 2 subsystem entanglement. 

In further considerations we shall show that thanks to the proper choice of the 
external bath parameters we are able to influence both: degree and character of the 
entanglement decay, leading to various interesting phenomena. Thus, in the present 
paper we shall concentrate on the squeezed vacuum state reservoir and on the influence 
of the squeezing parameters on the entanglement decay processes. 

To describe the evolution of the system influenced by a broadband squeezed vacuum 
we shall apply master equation approach for the reduced density matrix of the system. 
Including the standard rotating wave and Markov approximations we derive the master 
equation in the following form, analogously as in [lOj : 

tl = i(pH- Hp) + £ [2C jP C] - C)C jP - pC)C 3 



j=a,b 



+ 2 ]T [df^pCf + Cf ] pCf )] - Cf )] Cf ] p - pCfcf )] ] (4) 



=a,b 

E [A 

j=a.b 



Mt^t _ M^C^p - pCfc^ 



E [26 



j=a,b 



{s) pcf -cfcfp-p-cfcf 



where H = H n i + Hi + H e , the operators Cj, Cf^ and (j = a,b) describe the 
damping in the modes a and b, respectively, and are defined as: 

Cj = y/TjJ , (5) 
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= VT^ 3 > (6) 

Cf = y f^jj , (7) 

Mj = yjNjiNj + 1) exp(-i<t>) (8) 

The parameters jj are the spontaneous emission rates for both of the modes (j = {a, b}) 
, Mj characterises degree of squeezing of external bath, whereas <fi is the phase of 



squeezing for the bath field. If we assume that Mj = we deal with usual thermal 
bath characterized by the mean number of photons equal to Nj in the modes j = {a, b}. 
However, if we look at the form of eq.(8) the strength of the squeezing can be 
characterised by the value of Nj as well (and we shall use it in further considerations). 

As already pointed, the system populated under the interactions with external 
environment is more dimensional that considered at the beginning of the evolution 
(3 <g) 2 system). Anyway, we are interested in the dynamics within the subsystem 
which was occupied before the interactions were included. Therefore, to discuss the 
entanglement decay we extract 3(g>2 subspace occupied initially from the whole two-mode 
subspace which is involved in the system's dynamics, using the truncation procedure. 
Of course, one should keep in mind that entanglement can be temporally transferred 
to other subspace, different from the subspace defined by the states (eqnsHH). The 
density matrix for the truncated subspace can be obtained from the full density matrix 
calculated for 6x6. The truncation procedure is given by the following relation: 

Ptrunc = (|0)(0| + |1)(1| + |2)(2|)p(|0)(0| + |2)(2|) . (9) 

In this way we trace the evolution of entanglement which, before the interactions with 
environment, lived exclusively in the 3 <g) 2 space, keeping in mind that some amount 
of entanglement can be (and nodoubtly is) transferred periodically out of this traced 
subspace. 

Since our system is not a simple qubit-qubit one (for such a model the concurrence 
parameter for the entanglement description can be applied) but qutrit-qubit is 
considered, for determination of the degree of entanglement in the 3 <S> 2 subspace we 
use the negativity parameter j22j [23J [21]. It is defined using the eigenvalues (Xj) of the 
partial transpose, with respect to the smaller-dimensional subsystem, of the truncated 
system's density matrix (pj runc ) as follows: 

N{p) = max (0, -2A mi „) , (10) 

where \ m in is the smallest of the eigenvalues of the partially transposed density matrix 
of the system described. Negativity was proved to be an entanglement monotone [21] 
and as such it can be used as a good measure of entanglement. It takes a value of 1 for 
the MES states (for the case discussed here - Bell-like states) and for separable ones 
it is equal to zero. Using such an entanglement measure we are able to determine the 
degree of entanglement in the whole system considered. 

As already pointed out in [TU |25] , the system under consideration can be a source 
of maximally entangled states \B\), |-B 2 ) and \B$) {see eqns.flU - El)). It should be 



Squeezed vacuum reservoir effect for entanglement decay in nonlinear quantum scissors system.Q 

mentioned that although in the paper the entanglement dynamics in thermal and 
zero-temperature reservoirs has already been presented, however for the entanglement 
within two-qubit subspace, in contrast to the situation dealt with here, where the 3 <8> 2 
system is considered. As shown in [25], the entanglement sudden death and rebirth 
phenomena within the 2 <g> 2 subsystems considered were caused by the interaction of 
this subsystem with the third entangled state \B 3 ) from which the entanglement can be 
periodically transferred to the two qubit subsystem (formed by \Bi) and \B 2 ) states). 
We should point out that in the present paper we deal with entanglement generated in 
the qutrit-qubit (3 ® 2) system and we shall concentrate on the influence of squeezed 
reservoir on the entanglement decay process. 

For the system considered here we assume that the MES (Bell-like state \B%)) is 
generated initially and next, the external squeezed bath starts to influence the system. 
Since we are interested in the character of negativity decay when the system interacts 
with squeezed reservoir and especially, we want to check whether the squeezed reservoir 
can influence this decay in any way by examining the time-evolution of the negativity 
for various values of the parameters describing the reservoir squeezed field. As shown 
in our previous paper |25j, if the Bell state interacts with zero-temperature reservoir 
and there is no interaction between the system's parts forming entangled state, the 
entanglement decays asymptotically i.e. for t — > oo N(p) — > 0. However, if we deal 
with the thermal reservoir, the sudden death of entanglement phenomenon occurs in 
the system (entanglement decays to zero in a finite time). For the system interacting 
with the broadband squeezed vacuum bath (the case discussed here) we see that the 
moment of time of the last sudden death of entanglement phenomenon occurrence can 
depend on the squeezing parameters. In particular, for the system discussed here we 
assume that the bath is squeezed in the same mode as that in which the coupler is 
externally pumped. At first, looking at Fig. 2 we can see the influence of the squeezing 
parameter N a on final disentanglement time Td for constant squeezed vacuum phase 0. 
Similarly as for the case when we deal with not squeezed and zero-temperature bath 
(N a = 0) the time Td reaches high values - we can conclude that we obtain almost 
an asymptotic decay of entanglement. However, higher values of N a make the sudden 
death of entanglement visible as finite values of time Td become considerably smaller. 
Shortening of Td can be observed for both cases considered here (0 = {0, 7r}). Moreover, 
for these two cases the value of Td decreases continuously and tends asymptotically to 
its final value (see Fig. 2). However, comparing the two dependences of Td on N a , we 
can easily identify some values of N a for which is considerably different for = and 

= 7T. 

It means that by changing the value of we may prevent the system from 
rebuilding the entanglement despite the fact that there is still an interaction between 
both oscillators and a (strength of the interaction with external pumping field) is large 
enough to produce another entanglement birth. 

For further considerations we will choose non-zero values for N a and shall 
concentrate on the influence of the phase on the entanglement decay dynamics. In 
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particular, we shall find such a value of N a for which we can observe the noticeable 
changes in the values of as <fi varies form to 7T, concentrating in further considerations 
on the influence of the squeezed reservoir properties on rebuilding of the entanglement 
in the system. For instance, if we neglect mutual interaction between two oscillators 
(e = 0), assume weak external pumping (a = 0.02) and N a = 2, Nf, = (squeezed 
bath in one mode) sudden death of entanglement appears. We see (Fig. 3) that when 
the phase of squeezing is equal to zero (or it) the entanglement death appears slightly 
later than for the non-squeezed, thermal field case (M a = 0, N a = 2). However, if we 
additionally assume that apart from the absence of coupling between the oscillators the 
external coherent field is also cut off, the changes are practically unnoticeable and the 
result is practically identical to that of thermal reservoir case. Moreover, it is worth 
mentioning that any larger changes of the moment of time when the entanglement death 
appears are noticeable when the squeezing is present in the same mode of the reservoir 
modes that is externally pumped. 

At this point we assume that two oscillators interact with each other (i. e. e ^ 0) and 
we will focus on the differences in the entanglement decay due to the interactions with 
broadband squeezed vacuum, in comparison to the standard thermal reservoir. As shown 
in [25] the oscillations of entanglement appear in the system and for the thermal reservoir 
case besides the regular oscillations (characteristic for zero-temperature environment) 
it is possible to observe sequences of the sudden deaths and entanglement revivals. But 
what is more interesting, for the case of the squeezed external bath the phase of the 
squeezed reservoir mode can influence the time after which the final disentanglement is 
obtained and also influences the amount of the entanglement rebuilt in the system. In 
particular, we have noticed that by changing the phase of the squeezed reservoir mode 
we can qualitatively change the decay of negativity. It is seen for the whole range of e/a 
(for all cases: e/a > 1, e/a < 1 and e/a = 1). In particular, in Fig.4(a-c) we plotted the 
time after which the final entanglement death occurs (r^) as a function of squeezed bath 
phase showing that one additional rebuilt maximum in negativity versus time appears 
for some range of squeezed reservoir phases (~ 0; ~ 127r/20) and (~ 3l7r/20; ~ 2ir) (a 
similar feature can be observed for the thermal reservoir case, although there [25] the 
mean number of photons and external coupling strength dependences were discussed 
instead of phase dependence). For the system discussed here we can observe a significant 
shortening (~ 20%) of total disentanglement time for e/a = 1 (Fig. 4b) and shortening 
by ~ 17% for both e/a > 1 and e/a < 1 (see Figs. 4a and 4c). In Fig.2 we can see that 
this shortening (when we compare the cases for <fi = tt and = 0) varies from ~ 4% to 
~ 24% for various values of squeezing parameter N a . Our considerations concern one of 
these N a values that cause the maximal shortening. 

Moreover, it is seen that for all cases shown in Fig.4a-c for some values of <fi we 
observe one entanglement rebirth event less as a result of the influence of the squeezing 
phase on the entanglement dynamics. This change is of sudden character and indicates 
that we have some sort of the "phase transition" within the system discussed. For 
e/a > 1 and e/a < 1 (Fig.4a and Fig.4c) we can see that the interaction of the system 
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with the squeezed vacuum (for (f> G (~ 127r/20;~ 3l7r/20)) can shorten the time r d 
as we compare it with the case when ordinary thermal bath is considered. However 
this shortening is less pronounced - see the inset. Beyond this range even though a is 
not strong enough to produce another entanglement rebirth, the interactions with the 
squeezed bath lead to the appearance of additional negativity re-occurrence and hence, 
to a significant (~ 17%) lengthen of r d . When e/a — 1 (Fig. 4b) the changes of (ft are 
of the "phase transition" character again (as for the previous cases), but for this case 
one can observe vanishing of the last negativity maximum and consequently, a rapid 
shortening of T d time. 

The character of this "transition" is related to the fact that if the squeezing phase 
changes its value the terms in the master equation proportional to the exp(±i</>) add 
or subtract each other. In consequence, the populations and coherences in the system 
considered (and corresponding to them matrix elements of p) will be affected as well. 
Due to the fact that the definition of the negativity involves max function, changes 
in the matrix elements of p lead to the rapid changes (zeroing) in the negativity. A 
similar behaviour was discussed in [26] where the squeezing phase dependence of the 
decoherence rate was considered for the atomic system located inside the cavity. 

Thus, we can see that the phase of squeezed reservoir can influence not only the 
time of entanglement death but also can change the system's ability to rebuild the 
entanglement. Especially, we can see that even though the value of a is enough to 
cause the entanglement rebuild, the appropriately chosen phase <fi may prevent the last 
entanglement reconstruction. 

One more aspect that can be influenced by the squeezed reservoir state phase is 
the amount of entanglement rebuild in the system during the last entanglement rebirth 
phenomenon observed. Fig. 5 shows the dependence of the maximal value of the last 
M maXlast and last but one negativity maximum M m ax on the phase 0. Similarly, to 
the case discussed earlier (Fig. 4) we can observe some kind of "phase transition" - like 
behaviour - sudden changes of the maximal value of the last negativity reconstruction 
Af m ax la3t for some values of (ft (inset in Fig. 5) when the amplitude of the last negativity 
maximum is equal to zero. One should note that for the same range of </> we observe 
significant shortening of the time t& as well. Additionally, we can see from Fig. 5 that the 
interaction with the squeezed bath can reduce the amount of entanglement preserved in 
the system. This reduction depends on the value ofe/a>l-fore/a>l and e/a < 1 
it is about 30%, while for e/a = 1 even 56%. 

For the case of smaller squeezing parameter (iV a = 1, mode b is not squeezed) we 
can only observe the regular changes in the time after which entanglement disappears 
and we do not see any "phase transitions" in the system (Fig. 6). The shortening of the 
death time is only by about 4% and all the rebuild maxima of the negativity - time 
dependence are present for all squeezed reservoir phases — see Fig. 6. It confirms the 
earlier conclusion about the r d dependence on the squeezing parameter N a visible at 
Fig.2. It should be also mentioned that for this case the changes in value of the last 
rebuilt maximum are of the same character as for the case of stronger squeezing, but 
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smaller in percent. 

As we look at the time-evolution of both: negativity and populations of individual 
states, the question arises if the reservoir effects and population transfer could influence 
considerably the negativity evolution. The answer can be obtained for instance, from 
analysis of the trace of the density matrix ptmnc corresponding to the 3 <E> 2 subsystem 
considered. Thus, Fig. 7 shows the plots of the negativities and the trace of ptmnc for 
= and 7T. We see that for both values of </> the whole system that was initially 
populated in this subspace (Tr(p irimc ) = 1 for t=0) starts to populate other states - 
Tr(p trunc ) decreases and is practically identical for both values of <fi. The character 
of this decay depends on the value of the negativity for a given moment of time. In 
particular, the decay rate increases when the negativity goes to zero, whereas it slows 
down as the negativity reaches its maximum (this effect is visible in particular for the 
first negativity reconstruction). However those effects do not play crucial role, especially 
for longer times. This means that the population leakage can effect the entanglement 
dynamics to some extend. 

3. Summary 

We have analysed the influence of squeezed reservoir on the entanglement decay for a 
nonlinear system initially prepared in a maximally entangled state (one of the Bell-like 
state). Such a state determines some subspace of the states that can be treated as 
qutrit-qubit system and we have concentrated on it. Obviously, when the interaction 
with an external reservoir is included, other states are populated and in that way the 
whole system is no more a simple qutrit-qubit one, but more dimensional qudits should 
been considered. Anyway, our aim was to trace the evolution of the entanglement within 
the initially occupied subspace. Therefore, we look for the entanglement evolution inside 
the same qutrit-qubit system, keeping in mind that the other states (corresponding to 
the lower energies of the system) are involved in the dynamics and some loses of the 
entanglement in the subspace considered are related to the populating of these states. 
For the case when the oscillators could not interact with each other we could observe only 
small variations in the values of time when the entanglement disappears completely as 
a function of the squeezing phase. Moreover, this value does not differs considerably from 
that for the model involving standard thermal reservoir. However, the more pronounced 
changes become visible when both the interactions (external characterised by a and 
internal characterised by e) are present. From the discussion of a similar model [25] 
(where the thermal bath was considered) we expect oscillations of negativity, but their 
character depends on the squeezed reservoir properties. In particular, we have shown 
that shortening of the death time for the parameters used, by about 17% to 20% of its 
original value can be obtained while e/a changes for some range of squeezed reservoir 
phases. The changes in this shortening can be of the step character resembling in some 
sense a "phase transition" in the system. In fact, we observe one entanglement revival 
less or more, as a result of the interaction with the squeezed bath depending on the phase 
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of squeezing. Moreover, some changes in the amplitude of the last negativity maximum 
due to the squeezing phase variations could be also observed. In consequence, we can 
reduce the value of the entanglement revival at the end of its life by about 30% to 56% 
of its original value. Concluding, we can say that the squeezing parameters of external 
bath can influence considerably the entanglement decay dynamics leading to various 
interesting features. 
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Figure 1. The fidelity F between the two-mode state Y&) = c n,m{t) \n) a \m) b 

n.m—O 

and a state corresponding to the qutrit-qubit system l^)^ = Co2|0) |2)(,+ci2|l)o|2)b + 
C2o|2) a |0)b. The parameters: Xa = Xb = 25, a/e = 1, e = 0.1, -f a = lb = 0; time is 
scaled in 1/x units. 



350 




1 1 1 1 1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

N 

a 

Figure 2. Time of total entanglement death Td (scaled in 1/x units) versus the 
squeezing strength characterised by N a . Xa = Xb = 25, a — 0.01, e = 0.1, 
7a = 76 = 0.0025. We assume that = 0)) = \B 3 ) 
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Figure 3. Negativity N(p) versus time scaled in 1/x units; \ a — \b — 25, a = 0.01, 
e = 0, 7a = 7fc = 0.0025 and N a = 2, N b = 0, |*)jn = |B 3 >- 
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Figure 4. The dependence of as a function of the squeezing phase 4>. The system 
is initially prepared in a \B 3 ) state. The time is scaled in 1/x units, Xa = Xb = 25, 
e = 0.1 , 7„ = 7b = 0.0025 and N a = 2,N b = 0. At (a) e/a = 10, (b) ->• e/a = 1 and at 
(c) e/a = 0.5. Inset in (c) shows the central part of the whole plot. Stars correspond 
to thermal non-squeezed reservoir characterized by the same mean number of photons. 
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Figure 5. The maximal value of the last but one negativity reconstruction N max f° r 
the entanglement reborn appearance as a function of phase of squeezing reservoir </>. 
The system parameters are the same as in Fig. 4b. At the inset the amplitude of the 
last negativity maximum Af ma x last ■ 
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Figure 6. The same as in Fig. 4c but for iV a = 1 and JV& = 0. 
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Figure 7. Negativity for the qutrit-qubit subsystem described by the density matrix 
Ptrunc (eq.(9)) for </> = (solid line), <j) = it (dashed line) and the trace of ptrunc 
(dashed-dotted line). The remaining parameters are the same as in Fig. 4b. The trace 
of ptrunc is identical for both values of cf>. 



